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We compute the electromagnetic fields generated in heavy-ion collisions by using the HIJING 
model. Although after averaging over many events only the magnetic field perpendicular to the reac- 
tion plane is sizable, we find very strong electric and magnetic fields both parallel and perpendicular 
to the reaction plane on the event-by-event basis. We study the time evolution and the spatial distri- 
bution of these fields. Especially, the electromagnetic response of the quark-gluon plasma can give 
non-trivial evolution of the electromagnetic fields. The implications of the strong electromagnetic 
fields on the hadronic observables are also discussed. 

PACS numbers: 25.75.-q, 25.75.Ag 

I. INTRODUCTION 

Relativistic heavy-ion collisions provide us the methods to create and explore strongly interacting matter 
at high energy densities where the deconfined quark-gluon plasma (QGP) is expected to form. The prop- 
erties of matter governed by quantum chromodynamics (QCD) have been studied at the Relativistic Heavy 
Ion Collider (RHIC) at Brookhaven National Laboratory (BNL) and at the Large Hadron Collider (LHC) at 
CERN. Measurements performed at RHIC in Au + Au collisions at center-of-mass energy y/s = 200 GeV 
per nucleon pair and at LHC in Pb + Pb collisions at center-of-mass energy \/s = 2.76 TeV per nucleon pair 
have revealed several unusual properties of this hot, dense, matter {e.g., its very low shear viscosity [1,2], 
and its high opacity for energetic jets [3-6]). 

Due to the fast, oppositely directed, motion of two colliding ions, off-central heavy-ion collisions can 
create strong transient magnetic fields [7]. As estimated by Kharzeev, McLerran, and Warringa [8], the 
magnetic fields generated in off-central Au + Au collision at RHIC can reach eB ~ m 2 ~ 10 18 G, which 
is 10 13 times larger than the strongest man-made steady magnetic field in the laboratory. The magnetic 
field generated at LHC energy can be 10 times larger than that at RHIC [9]. Thus, heavy-ion collisions 
provide a unique terrestrial environment to study QCD in strong magnetic fields. It has been shown that a 
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strong magnetic field can convert topological charge fluctuations in the QCD vacuum into global electric 
charge separation with respect to the reaction plane [8, 10, 1 1]. This so-called chiral magnetic effect may 
serve as a sign of the local P and CP violation of QCD. Experimentally, the STAR [12, 13], PHENIX [14], 
and ALICE [15] Collaborations have reported the measurements of the two-particle correlations of charged 
particles with respect to the reaction plane, which are qualitatively consistent with the chiral magnetic effect, 
although there are still some debates [16-20]. 

Besides the chiral magnetic effect, there can be other effects caused by the strong magnetic fields in- 
cluding the catalysis of chiral symmetry breaking [21], the possible splitting of chiral and deconfinement 
phase transitions [22], the spontaneous electromagnetic superconductivity of QCD vacuum [23, 24], the 
possible enhancement of elliptic flow of charged particles [25, 26], the energy loss due to the synchrotron 
radiation of quarks [27], the emergence of anisotropic viscosities [26, 28, 29], the induction of the electric 
quadrupole moment of the QGP [30], etc. 

The key quantity of all these effects are the strength of the magnetic fields. Most of the previous works 
estimated the magnetic-field strength based on the averaging over many events [8, 9, 31-33], thus, due to 
the mirror symmetry of the collision geometry, only the y-component of the magnetic field remains sizable, 
e(B y ) ~ m^, while other components are (B x ) = (B z ) = 0. Hereafter, we use the angle bracket to denote 
event average. We choose the z axis along the beam direction of the projectile, x axis along the impact 
parameter b from the target to the projectile, and y axis perpendicular to the reaction plane, as illustrated in 
Fig. 1. 
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FIG. 1: The geometrical illustration of the off-central collisions with impact parameter b. Here "T" for target and "P" 
for projectile. 

However, in many cases, the final hadronic signals are measured on the event-by-event basis. Thus, it 
is important to study how the magnetic fields are generated on the event-by-event basis. Such a study was 
recently initiated by Bzdak and Skokov [34]. They showed that the x-component of the magnetic field can 
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be as strong as the y component on the event-by-event basis, due to the fluctuation of the proton positions 
in the colliding nuclei. Besides, they also found that the event-by-event generated electric field can be 
comparable to the magnetic field. 

The aim of our work is to give a detailed study of the space-time structure of the event-by-event generated 
electromagnetic fields in the heavy-ion collisions. We perform our calculation by using the heavy ion 
jet interaction generator (HIJING) model [35-37]. HIJING is a Monte-Carlo event generator for hadron 
productions in high energy p + p, p + A, and A + A collisions. It is essentially a two-component model, 
which describes the production of hard parton jets and the soft interaction between nucleon remnants. In 
HIJING, the hard jets production is controlled by perturbative QCD, and the interaction of nucleon remnants 
via soft gluon exchanges is described by the string model [38]. 

This paper is organized as follows. We give a general setup of our calculation in Sec. II. In Sec. Ill, 
we present our numerical results. We discuss the influence of the non-trivial electromagnetic response of 
the QGP on the time evolution of the electromagnetic fields in Sec. IV. We conclude with discussions and 
summary in Sec. V. We use natural unit h = c = 1. 

II. GENERAL SETUP 

We use the Lienard-Wiechert potentials to calculate the electric and magnetic fields at a position r and 
time t, 



where Z n is the charge number of the nth particle, R n = r — r n is the relative position of the field point 
r to the source point r n , and r n is the location of the nth particle with velocity v n at the retarded time 
t n = t — |r — r n |. The summations run over all charged particles in the system. Although there are 
singularities at R n = in Eqs. (2.1)-(2.2), in practical calculation of E and B at given (i,r), the events 
causing such singularities rarely appear. So, we omit such events in our numerical code. In non-relativistic 
limit, v n <C 1, Eq. (2.1) reduces to the Coulomb's law and Eq. (2.2) reduces to the Biot-Savart law for a set 
of moving charges, 




(2.1) 



(2.2) 




(2.3) 



(2.4) 
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To calculate the electromagnetic fields at moment t, we need to know the full phase space information 
of all charged particles before t. In the HIJING model, the position of each nucleon before collision is 
sampled according to the Woods-Saxon distribution. The energy for each nucleon is set to be y/s/2 in the 
center-of-mass frame. Assuming that the nucleons have no transverse momenta before collision, the value 
of the velocity of each nucleon is given by v z = 1 — (2rriN / '\fs) 2 , where is the mass of the nucleon. 
At RHIC and LHC, v z is very large, so the nuclei are Lorentz contracted to pancake shapes. 

We set the initial time t = as the moment when the two nuclei completely overlap. The collision time 
for each nucleon is given according to its initial longitudinal position zjy = z^ ■ 2m^ / 'yfs and velocity v z , 
where is the initial longitudinal position of nucleon in the rest frame of the nucleus. The probability 
of two nucleon colliding at a given impact parameter b is determined by the Glauber model. In this paper, 
we call nucleons without any interaction spectators and those that suffer at least once elastic or inelastic 
collision participants before their first collision. The participants will exchange their momenta and energies, 
and become remnants after collision. Differently from the spectators and the participants, the remnants can 
have finite transverse momenta. In our calculation, we find that although the spectators and participants are 
the main sources of fields at t < 0, remnants can give important contributions at t > 0. In the HIJING 
model, we neglect the back reaction of the electromagnetic field on the motions of the the charged particles. 
This is a good approximation before the collision happens because the electromagnetic field is weak at that 
time. We will discuss the feed back effect of the electromagnetic fields on QGP in Sec. IV by using a 
magnetohydrodynamic treatment. 

After collision, many partons may be produced and the hot, dense, QGP may form. As the QGP is nearly 
neutral, we neglect the contributions from the produced partons to the generation of the electromagnetic 
field. However, if the electric conductivity of the QGP is large, the QGP can have significant response to 
the change of the external electromagnetic field. This can become substantial for the time evolution of the 
fields in the QGP. We will discuss this point in detail in Sec. IV. 

III. NUMERICAL RESULTS 
A. Impact parameter dependence 

We first show the impact parameter dependence of the electromagnetic fields at r = and t = 0. The 
left panel of Fig. 2 is the results for Au + Au collision at RHIC energy \fs = 200 GeV; the right panel of 
Fig. 2 is for Pb + Pb collision at LHC energy y/s = 2.76 TeV As seen from Eq. (2.2), (B x (t, 0)) = 0, 
while (B v (t, 0)) < when b > 0. Also, from Eqs. (2.1)-(2.2), we find that there are always \E y (0, 0)| « 
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FIG. 2: (Color online) The electromagnetic fields at t = and r = as functions of the impact parameter b. 

\B X (0,0)\ and \B y (0,0)\ > \E X (0,0)\ when v z is large [See Eqs. (3.2)-(3.3)]. These facts are reflected 
in Fig. 2. Although the x-component of the magnetic field as well as the x- and y-components of the 
electric field vanish after averaging over many events, their magnitudes in each event can be huge due to 
the fluctuations of the proton positions in the nuclei. Thus, following Bzdak and Skokov [34], we plot 
the averaged absolute values (\E Xjy \) and (\B x . y \) at r = and t = 0. Similar with the findings in 
Ref. [34], we find that {\B X \}, (\E X \), and {\E y \} are comparable to {\B y \}, and the following equalities 
hold approximately, ~ (l^yl) ~ (l-^U)- But our results at RHIC energy are about three times 

smaller than that obtained in Ref. [34]. We checked that this is because the thickness of the nuclei in our 
calculation is finite while the authors of Ref. [34] assumed that the nuclei are infinitely thin. We can also 
observe that, at small b region, contrary to (B y ) which is proportional to b, the fields caused by fluctuation 
are not sensitive to b. 



We see from Fig. 2 that the magnitudes of all the fields at LHC energy is around 14 times bigger than 
that at RHIC energy. To study the collision energy dependence more carefully, we calculate the fields at 
t = and r = for different yfs. To high precision, the linear dependence of the fields on the collision 
energy is obtained, as shown in Fig. 3. Thus, the following scaling law holds for event-by-event generated 
electromagnetic fields as well as for event-averaged magnetic fields, 



where Ra is the radius of the nucleus and is a universal function which has the shapes as shown 



B. Collision energy dependence 



e • Field oc y/sf(b/R A ) 



(3.1) 



in Fig. 2 for ( | B x>y \),{\E x>y \), and (B y ). 
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FIG. 3: (Color online) The collision energy dependence of the electromagnetic fields at r = and t = 0. 



Actually, a more general form of Eq. (3.1) can be derived from Eqs. (2.1)-(2.2). As the fields at t = 
are mainly caused by spectators and participants whose velocity v n = v z = \Jl — (2mjv/y / s) 2 « 1, the 
electric and magnetic fields at t = in the transverse plane can be expressed as 

^(0,r) « f ^£^, ( 3 - 2 ) 
4vr 2m^ ^ R 3 nl 

eB ± (0,r) « — > -3 , (3.3) 

4vr 2m w ^ R* ± 

where e nz is the unit vector in ±z direction depending on whether the nth proton is in the target or in the 
projectile, R n j_ is the transverse position of the nth proton which is independent of y/s, and R n ± = |R n j_|. 

C. Spatial distribution 

The spatial distributions of the magnetic and electric fields are evidently inhomogeneous. We show in 
Fig. 4 the contour plots of (B Xj y :Z ), (E Xj y jZ ), (\B X;y>z \), and (\E Xj y :Z \) at t = in the transverse plane 
at RHIC energy. The upper two panels are for b = and the lower two panels are for b = 10 fm. The 
spatial distribution of the transverse fields for LHC energy is merely the same as Fig. 4 but the fields have 
2760/200 as 14 times larger magnitudes everywhere according to Eqs. (3.2)-(3.3). The spatial distribution 
of the fields in the reaction plane was studied in Ref. [32]. 
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FIG. 4: (Color online) The spatial distributions of the electromagnetic fields in the transverse plane at t = for b = 
(upper panels) and b — 10 fm (lower panels) at RHIC energy. The unit is mj. The dashed circles indicate the two 
colliding nuclei. 



First, as we expected, the longitudinal fields (B z ), (E z ), (\B Z \), and (\E Z \) are much smaller than 
the transverse fields. Second, the event-averaged fields (B XjV ) and (E xy ) distribute similarly with the 
fields generated by two uniformly charged, oppositely moving, discs. Third, the spatial distribution of the 
magnetic fields is very different from that of the electric fields on the event-by-event basis. For central 
collisions, both (l-B^I) and {\B y \} distribute circularly and concentrate at r = 0, while (l-E^I) and (\E y \) 
peak around x = ±Ra and y = ±Ra with Ra the radius of the nucleus. We notice that for off-central 
collision, the y-component of the electric field varies steeply along y-direction, reflecting the fact that at 
t = a large amount of net charge stays temporally in the "almond"-shaped overlapping region. 



D. Probability distribution over events 

Although we used the event-averaged absolute values (\B Xjy \) and (\E x . y \) to characterize the event-by- 
event fluctuations of the electromagnetic fields, it would have more practical relevance to see the probability 
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FIG. 5: (Color online) The probability densities P(B X , B y ) and P(E x ,E y ) for different impact parameters for Au + 
Au collisions at y/a = 200 GeV. 



P{B x ,B y ) = - — — , (3.4) 



distribution of the magnetic field, denned as 

1 d 2 N 
N dB x dB y 

where N is the number of events. Similarly, we can define P(E X , E y ). After simulating 10 6 events, we 
obtain P(B X , B y ) and P(E X , E y ) at t = and r = for Au + Au collisions at ^/s = 200 GeV, as shown in 
Fig. 5. As expected, the probability distribution of the magnetic (electric) field peaks at B = (E = 0) for 
central collision, while the probability distribution for magnetic field is shifted to finite B y for off-central 
collision. This is more clearly shown in Fig. 6, where we depict the one-dimensional probability density 
P(B X ) = f dB y P(B x , B y ) (other probability densities are analogously defined). 

The probability distributions for Pb + Pb collisions at = 2.76 TeV have analogous shapes with Fig. 5 
but much more spread, as clearly shown in the lower panels of Fig. 6. This is because the strength of the field 
generated at LHC can be obtained approximately from that at RHIC by a \J slhc / •SRHic- scaim g according 
to Eqs. (3.2)-(3.3). Thus, after normalization, the probability distributions at LHC energy are related to that 
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at RHIC energy by 
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FIG. 6: (Color online) The probability densities P(B X:V ) and P{E XtV ) for central collision 6 = and off-central 
collision b = 10 fm for Au + Au collision (upper panels) at s/s = 200 GeV and for Pb + Pb collision at s/s = 2.76 
TeV (lower panels). 



E. Early-stage time evolution 

In Fig. 7, we show our results of the early-stage time evolution of the electromagnetic fields at r = 
in both central collisions and off-central collisions with b = 10 fm, for Au + Au collision at y/s = 200 
GeV and for Pb + Pb collision at yfs = 2.76 TeV. We take into account the contributions from charged 
particles in spectators, participants, and remnants. Around t = 0, we checked that the contributions from 
the remnants are negligibly small, while the contributions from participants can be as large as that from 
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FIG. 7: (Color online) The time evolution of electromagnetic fields at r = with impact parameter b = 0, and 6=10 
for Au + Au collisions at s/s = 200 GeV and Pb + Pb collisions at s/s = 2.76 TeV. After collision, the remnants can 
essentially slow down the decay of the transverse fields, and enhance the longitudinal fields. 



spectators. However, at a latter time when the spectators have already moved far away from the collision 
region, the contributions from the remnants become important because they move much slower than the 
spectators. These remnants can essentially slow down the decay of the transverse fields, as seen from Fig. 7. 
Another evident effect of the remnants is the substantial enhancements of the longitudinal magnetic and 
electric fields caused by the position fluctuation of the remnants, which have non-zero transverse momenta. 
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Particularly, although (|-B Z |) is at least one order smaller than (\B Xty \) for t < 1 fm/c, (\E Z \) can evolve to 
the same amount of (l-E^I) in a very short time after collision and then they decay very slowly. 

For central collision, all the fields are generated due to the position fluctuations of the charged particles. 
As seen from Fig. 7, these fluctuations lead to sizable (|£? x |) = {\E y \) and (|-B 2 |) = (\B y \) around t = 0, 
but they drop very fast. Note that the fields drop faster for larger collision energy yfs. 

For off-central collisions, the y-component of the magnetic field are much larger than other fields at 
t = 0. But at a latter time when the spectators move far away, the contributions of remnants dominate, and 
leadto (1^1) « (1^1) 

A common feature of all the fluctuation-caused transverse fields is that they all increase very fast before 
collision (due to the fast approaching of the nuclei), then they drop steeply after t = (due to the high- 
speed leaving of the spectators away from the collision center), and then decay very slowly (due to that 
contribution from the slowly moving remnants take over that from the spectators). After the early-stage 
evolution, the produced QGP may get enough time to respond to the electromagnetic fields, which may 
substantially modify the picture of evolution. We discuss this point in the next section. 



IV. RESPONSE OF THE QGP TO ELECTROMAGNETIC FIELDS 

In the calculations above, we have neglected the electromagnetic response of the matter produced in the 
collision (i.e., we assumed the produced matter is ideally electrically insulating). However, if the produced 
matter, after a short early-stage evolution, is in the QGP phase, the electric conductivity a is not negligible. 
At high temperature, the perturbative QCD gives that a ss 6T/e 2 [39], and the lattice calculations give 
that a ps 7C E mT [40], or a « OACemT [41, 42], or a sa (1/3)C EM T-C EM T [42, 43] for temperature 
of several T c , where Cem = Ylf e /> / = u i d, s > an( ^ e f * s tne charge of quark with flavor /. Thus it is 
expected that the QGP can have non-trivial electromagnetic response. Such electromagnetic response can 
substantially influence the time evolution of the electromagnetic fields in the QGP. 

To have an estimation of the electromagnetic response of QGP, we use the following Maxwell's equa- 
tions, 

dB 

VxE = - — , (4.1) 
ot 

-VxB = e^+J, (4.2) 
H at 

where fi and e are the permeability and permittivity of the QGP, respectively, and are assumed as constants. 
J is the electric current determined by the Ohm's law, 

J = a (E + v x B) , (4.3) 



12 



where v is the flow velocity of QGP. Using Eq. (4.3), we can rewrite the Maxwell's equations as 

f =Vx(«B) + I(VB-^), 



« ■ 8( xB = vx(VxE) + -^E-^, (4.5) 

where we have used the Gauss's laws V • B = and V • E = p = with the assumption that the net electric 
charge density of the QGP is zero. Equation (4.4) is the induction equation, which plays a central role in 
describing the dynamo mechanism of stellar magnetic field generation. The first terms on the right-hand 
sides of Eqs. (4.4)-(4.5) are the convection terms, while the last terms are the diffusion terms. The ratio of 
these two types of terms are characterized by the magnetic Reynolds number R m , 

R m = LUap, (4.6) 

where L is the characteristic length or time scale of the QGP, U is the characteristic velocity of the flow. 

Because the theoretical result of a is quite uncertain, the value of R m is also uncertain. For example, by 
assuming p = e = 1, setting the characteristic length scale L ~ 10 fm, and the typical velocity U ~ 0.5, 
we can estimate R m atT = 350 MeV as R m ~ 0.2 if we use a « OACemT in Refs. [41, 42], or R m ~ 4 
if we use a « 7C E mT in Ref. [40], or R m ~ 600 if we use a « 6T/e 2 in Ref. [39]. 

If R m <C 1, we can neglect the convection terms in Eq. (4.4) and Eq. (4.5). Tuchin studied this case [27] 
(with additional condition U <C 1 so that the second-order time-derivative terms in the diffusion terms are 
neglected), and concluded that the magnetic field can be considered as approximately stationary during the 
QGP lifetime 1 . 

If R m S> 1, we can neglect the diffusion terms in Eqs. (4.4)-(4.5), i.e., take the ideally conducting limit, 

— = Vx(vxB), (4.7) 
E = -v x B. (4.8) 

It is well-known that Eq. (4.7) leads to the frozen-in theorem for ideally conducting plasma {i.e., the mag- 
netic lines are frozen in the plasma elements or more precisely the magnetic flux through a closed loop 
defined by plasma elements keeps constant [44]). 

We now use Eqs. (4.7)-(4.8) to estimate how the electromagnetic field evolves in a QGP with R m 3> 1. 
To this purpose, we have to know the evolution of v first. By assuming that the bulk evolution of the QGP 



1 To reach this conclusion, Tuchin used a — 6T 2 /T c to estimate the magnetic diffusion time r = (L/2) 2 a/4 and found that, 
for L = 10 fm and T — 2T C w 400MeV, r w 150 fm is much longer than the lifetime of the QGP. However, if, for example, 
cr « OACemT is used, the magnetic diffusion time is r ~ 0.3 fm, which is much shorter than what Tuchin obtained. 
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is governed by strong dynamics, we can neglect the influence of the electromagnetic field on the evolution 
of the velocity v. We assume the Bjorken picture for the longitudinal expansion, 

z 



t 



(4.9) 



Because the early transverse expansion is slow, following Ref. [45], we adopt a linearized ideal hydrody- 
namic equation to describe the transverse flow velocity vi , 



1 



:V_l 



-y±P = -cJV± Ins, 



(4.10) 



dt'^ e + P 

where we used e + P = Ts, s is the entropy density, and c s = \J dP/de is the speed of sound. For 
simplicity, we choose an initial Gaussian transverse entropy density profile as in [45], 



s(x, y) = s exp 



_a_ V_ 



(4.11) 



where a xy are the root-mean-square widths of the transverse distribution. They are of order of the nuclei 
radii if the impact parameter is not large. For example, for Au + Au collisions at RHIC, a x ~ a y ~ 3 fm 
for 6 = 0, and a x ~ 2 fm, a y ~ 3 fm for 6 = 10 fm. One can then easily solve Eq. (4. 10) and obtain, 



n "fit) 
1^ 



yt. 



(4.12) 
(4.13) 



Substituting the velocity fields above into Eq. (4.7), we obtain a linear differential equation for B(i). For a 
given initial condition B°(r) = B(t = to, r) where is the formation time of the QGP, it can be solved 
analytically, 



B x (t,x,y,z) = je 2a y 



, ye y , z— 



B y (t,x,y,z) 



to _-3r(* 2 -*§) D o 



B" y xe 2a >- 



-(t 2 -® -^(* 2 ~*o) to 
,ye y 



(4.14) 



(4.15) 



t M ' t I 

Because B z is always much smaller than B x and B y , we are not interested in it. The electric fields can be 
obtained from Eq. (4.8) once we have B(t, r). 

To reveal the physical content in Eqs. (4.14)-(4.15), we notice that, by integrating Eq. (4.9), Eq. (4.12), 
and Eq. (4. 13), a fluid cell located at (xq, yo, zq) at time to will flow to the coordinate (x, y, z) at time t with 



x = xo exp 

V = yo exp 
t 



Z = Zo 



to' 



(4.16) 
(4.17) 
(4.18) 
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Thus,we can rewrite Eqs. (4.14)-(4.15) as 



t -A^-tD 



B x (t,x,y,z) = —e 2a y B x (t ,x ,y , z ) , (4.19) 
to c \ (t 2 f2 ) 

B y (t,x,y,z) = —e 2a * B y (t , x , y , z ) . (4.20) 



As the areas of the cross sections of the QGP expand according to t exp \ in the y-z plane and 

t exp (j^zt 2 ^j in the x-z plane, Eqs. (4.19)-(4.20) mean that the magnetic line flows with the fluid cell and 
is diluted due to the expansion of the QGP. These are just the manifestations of the frozen-in theorem, We 
also note that Eqs. (4.14)-(4.15) can be written in explicit scaling forms, 

° 2 ^2 / c s 2 c s £2 \ c s £2 / c s 2 c s ^2 \ 

te 2 ^ B x it,xe^ ,ye 2a y ,zt\ = t e 2 ^ °B X it , xe^ , ye 2 ^ °, zt , (4.21) 

c s t 2 ( C ° t 2 -^rt 2 \ C *t 2 { C » t 2 -^rt 2 \ 

te 2 ^ Bylt.xe 2 ^ ,ye 2a y ,zt \ = t e 2 ^ °B y [to, xe 2 ^ °, ye 2a y , z t • (4.22) 



From Eqs. (4.14)-(4.15), we see that the evolution of B is strongly influenced by its initial spatial distri- 
bution. However, the time evolution of the magnetic fields at the center of the collision region, r = 0, takes 
very simple forms, 

B x (t,0) = %~^ {e ~ tl) Bl{Qi), (4.23) 
By (t,0) = ^e'^'^B^O). (4.24) 

Setting a x ~ a y ~ 3 fm and c 2 s ~ 1/3, we see from Eqs. (4.23)-(4.24) that for t < 5 fm the magnetic fields 
decay inversely proportional to time. 



V. SUMMARY AND DISCUSSIONS 



In summary, we have utilized the HIJING model to investigate the generation and evolution of the 
electromagnetic fields in heavy-ion collisions. The cases of Au + Au collision at -y/s = 200 GeV and Pb + 
Pb collision at yfs = 2.76 TeV are considered in detail. Although after averaging over many events only 
the component B y remains, the event-by-event fluctuation of the positions of charged particles can induce 
components B x , E x , E y as large as B y . They can reach the order of several m 2 /e. The spatial structure of 
the electromagnetic field is studied and a very inhomogeneous distribution is found. We study also the time 
evolution of the fields including the early-stage and the QGP-stage evolutions. We find that the remnants can 
give substantial contribution to the fields during the early-stage evolutions. The non-trivial electromagnetic 
response of the QGP, which is sensitive to the electric conductivity, gives non-trivial time dependence of 
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the fields in it, see Sec. IV. We check both in numerical calculation (Fig. 3) and analytical derivations 
(Eqs. (3.2)-(3.3)) that the electric and magnetic fields at t = have approximately linear dependence on the 
collision energy y/s. 

The strong event-by-event fluctuation of the electromagnetic field may lead to important implications 
for observables which are sensitive to the electromagnetic field. We point out two examples here. 

(1) From Fig. 2 we see that although the electric fields E x and E y at r = can be very strong, they 
are roughly equal. Then one expects that the strong electric fields should not have a significant contribution 
for the correlation observable (cos(0i + 4>2~ 2^pp)} which is sensitive to the chiral magnetic effect [46], 
where <p\^ are the azimuthal angles of the final-state charged particles, and ^rp is the azimuthal angle 
of the reaction plane. On the other hand, from Fig. 4 we see that in the overlapping region for peripheral 
collisions, the electric field perpendicular to the reaction plane has larger gradient than that parallel to the 
reaction plane. Thus, a strong, out-of-plane electric field develops away from the origin r = in the 
overlapping region. Note that the direction of this electric field points outside of the reaction plane. If 
the electric conductivity of the matter produced in the collision is large, this out-of-plane electric field can 
drive positive (negative) charges to move outward (toward) the reaction plane, and thus induce an electric 
quadrupole moment. Such an electric quadrupole moment, as argued in Ref. [30], can lead to an elliptic 
flow imbalance between ir + and 7r _ . Thus, it will be interesting to study how strong the electric quadrupole 
moment can be induced by this out-of-plane electric field. Note that such electric quardrupole configuration 
does not contribute to the correlation (cos(0i + 4>2 — 2^rp)). 

(2) It is known that quarks produced in off-central heavy-ion collision can be possibly polarized due 
to the spin-orbital coupling of QCD [47-49]. The strong magnetic field can cause significant polarization 
of quarks as well through the interaction between the quark magnetic moment and the magnetic field. As 
estimated by Tuchin [27], a magnetic field of order m^/e can almost immediately polarize light quarks. 
Such polarization, contrary to the polarization due to spin-orbital coupling, will depend on the charges 
of quarks, and build a spin-charge correlation for quarks, i.e., the positively charged quarks are polarized 
along the magnetic field while the negatively charged quarks are polarized opposite to the magnetic field. 
If we expect that the strong interaction in the QGP and the hadronization processes do not wash out this 
spin-charge correlations, we should observe similar spin-charge correlation for final-state charged hadrons. 
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